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We report intrinsic tunnelling data for mesa structures fabricated on three over- and optimally- 
doped Bi2.i5Sri.s5CaCu20g+6 crystals with transition temperatures of 86-80 K and 0.16-0.19 holes 
per CuOa unit, for a wide range of temperature (T) and applied magnetic field ( H ), primarily 
focusing on one over-doped crystal. The differential conductance above the gap edge shows clear 
structure with informative T and H dependence. Data below the gap edge imply that the tunnelling 
process is mainly coherent and give clear evidence for strong T-dependent pair breaking. These 
findings are important for a detailed analysis in terms of Eliashberg theory and could help to clarify 
the pairing mechanism. 


Ever since the discovery of cuprate superconductors 
over 25 years ago, understanding the pairing mechanism 
has been an extremely active area of research. For clas¬ 
sical strong-coupling superconductor such as Pb, precise 
measurement of the current-voltage (/ — V) character¬ 
istics of planar tunnel junctions with a normal metal, 
or between two superconductors (SIS), shows good cor¬ 
relation between structure in dl/dV or d 2 I/d,V 2 above 
the superconducting (s/c) energy gap A, and the energy- 
dependent phonon density of states (DOS) measured by 
neutron scattering. Analysis using Eliashberg theory 
shows clearly that pairing is mediated by the virtual ex¬ 
change of phonons [T]. The structure in dl / dV is caused 
by the energy (or frequency ui) dependence of two com¬ 
plex parameters: A(w), associated with s/c order and 
z a (ui), arising from many-body, electron-phonon renor¬ 
malization in the s/c state that are coupled via two in¬ 
tegral equations. In the d- wave case they will also have 
strong momentum (k) dependence. 

Mesa structures fabricated from highly anisotropic 
high-T c superconductors such as Bi 2 Sr 2 CaCu 2 08+<5 (Bi- 
2212) may be regarded as stacks of planar “intrinsic tun¬ 
nel junctions” (ITJs) connected in series, and their I — V 
characteristics correspond to c-axis, SIS tunnelling spec¬ 
tra [2|. Apart from a very recent paper |3j, ITJs have not 
been used to detect structure in dl/dV that could give 
information about possible boson energies [3] and there 
are several reasons for this, see Footnote E3 However, 
tunnelling in planar ITJs may be easier to understand 
than in break junctions 0 □ and being a bulk probe 
it could be a valuable complement to angle-resolved pho¬ 
toemission spectroscopy (ARPES) [Sj , optical reflectivity 
g] and scanning tunnelling microscopy (STM) fijl-fTTl . 

Here we report experimental data showing structure in 
the dl/dV curves of ITJs fabricated on two over-doped 
single crystals of Bi-2212 with T c values of 80 and 78 K, 
OD80 and OD78, and an optimally doped crystal, OP86 


with T c = 86 K. We focus on OD80, for which there is 
clear structure above 2A 0 , and where ARPES data gives 
evidence for a large Fermi surface and no pseudogap at 
low T [5]. Data for hole concentrations p < 0.19 per 
Cu0 2 unit, may be complicated by the presence of the 
pseudogap and also of charge density waves that have 
been observed recently for both Bi-2212 m nn and 

YBa 2 Cu306.7 [HI- 

Single crystals of Bi-2212 were grown using a travelling 
solvent floating zone furnace and feed rods with nominal 
stoichiometry of Bi 2 ^Sr! 8 5 CaCu 2 08 + <5. These have a 
maximum T c of 86.5 K and we infer p from the empirical 
relation T c = T™ ax (\ — 82.6[p — 0.16] 2 ) [T3|, finding p = 
0.194, 0.191 and 0.16 for the three crystals studied. For 
OD78 and OD80, T c measured by Squid magnetometry 
before fabrication of mesas agreed with values deduced 
from mesa dl/dV curves, e.g. Fig.[2](b), to within 1 K but 
for OP86, T c from dl/dV curves could be up to 3 K lower 
than the magnetometry data. Typical I — V character¬ 
istics for the three mesas at small bias, measured when 
increasing /, are given as supplemental material (SM). 
The branches in the I—V characteristic correspond to 
different numbers of Josephson junctions being switched 
into the resistive state. The number of junctions ( N ) 
in the stack is equal to the total number of branches ob¬ 
served. In the SM it is shown that these branches scale on 
to each other, confirming that the junctions in the mesa 
have uniform area, and therefore all junctions switched 
to the resistive state will have the same voltage bias. 

Mesa dl/dV spectra pH] were measured as the bias 
current was swept down from its maximum value towards 
zero, thereby maintaining the resistive state. A lock-in 
technique with a small 77.7 Hz current modulation was 
employed. Fig. [lja) shows dl/dV per unit area for a sin¬ 
gle junction of OD80, at 1.4 K vs. the bias voltage per 
junction, for H = 0 to 13T applied perpendicular to the 
Cu0 2 planes. The curves are symmetric for ±E, so for 
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FIG. 1: Color online: (a) dl/dV curves for OD80 at 1.4 K 
in various magnetic fields applied perpendicular to the CuC >2 
planes plotted vs. V/N, the bias voltage per junction, (b) and 
(c) show details of the field dependence of the lower and upper 
dips. Here dl/dV curves have been normalized by dividing 
through by the polynomial given in the text. Values of (Ao + 
)/e, (2Ao + S2)/e and 2(Ao + H)/e are shown by arrows (see 
text). 


clarity we only show data for V > 0. For such SIS junc¬ 
tions the sharp peaks are located at voltages of 2A 0 /e, 
where Ao is the maximum value of the d -wave gap. At 
1.4 K this gives A 0 = 27.3, 26.9 and 27.4 meV for OD78, 
OD80 and OP86 respectively, in good agreement with the 
lower values shown in Fig. 15 of Ref. [9] for these doping 
levels. The ratio 2Ao/A:sT c = 8.08 ±0.1, 7.83 ±0.1 and 
7.5±0.15 for these three mesas is ^1.75-1.9 times larger 
than for a weak-coupling d -wave superconductor. 

In order to see whether tunnelling is coherent, i.e. 
whether it conserves k parallel to the CuC >2 planes, we 
make a log-log plot of dl/dV - (dI/dV) res vs. V, where 
(dl/dV)res is the extrapolated value at V = 0. As shown 
in the SM, for all three mesas this gives a limiting V 4 
law for V < 0.74Ao, whereas a V 2 law is expected for 
completely incoherent tunnelling. It should be possible 
to account for this V 4 law by including both the angle- 
dependence of the c-axis transfer integral, which goes to 
zero in the nodal directions m , and the Dynes damping 
factor [16], i.e. the imaginary part of A(w). As men¬ 
tioned already, this is likely to be much more w and k 
dependent than in the s-wave case. Defects in the barrier 
regions could scatter the tunnelling electrons sufficiently 
to give limited incoherence, i.e. a spread in k values m 
and this would be most significant in the nodal regions. 

The data for OD80 in Figs. [lja)-(c) , show two field (H) 
dependent dips above eV = 2A 0 . For clarity of display, 
the data in Figs. 0b) and (c) have been normalized by 
dividing through by the polynomial (8.85 — 151.27V 2 ± 
501.22V 4 ) x 10where V is the voltage per 


junction, see Footnote H8T 

For coherent tunnelling, the angle 9 between k and 
the anti-nodal direction is conserved. Within the stan¬ 
dard “semiconductor” model [lj , for a d- wave energy gap 
varying as Ao cos(20) and a single k-independent boson 
energy D, boson-induced structure is expected to be most 
apparent at eV=2A 0 ± S2. At this bias voltage, the gap 
edge at A 0 for 9 = 0 on one side of the junction has the 
same energy as any structure at A u ± O and 9 = 0 on the 
other side. The effect is larger there because the peak in 
the s/c quasi-particle DOS is largest at 9 = 0, with fur¬ 
ther enhancement from the larger normal state DOS [19] 
and tunnelling matrix elements [15] , Additional struc¬ 
ture is expected near eV=2A 0 ± 20 where boson-induced 
anomalies on each side of the junction at 9 = 0 have the 
same energy. The 5=1, magnetic resonance excitation, 
seen by inelastic neutron scattering [20], is a candidate 
pairing boson m E3. It has an energy 0 = 5.41 cbT c 
QHH] and a momentum vector Q, close to (n/a,Tr/a), 
na uni where a is the in-plane lattice spacing. Vari¬ 
ous energies associated with this value of O are shown 
in Fig. [I] and later for all three mesas in Fig. [3] It can 
be seen that there is only rough correspondence with the 
simple description given above. There is also strong H 
and T-dependence as amplified below. 

The quasi-particle DOS produced by a magnetic field 
in a d- wave superconductor, is predicted j22j to be ~ 
N( 0)(U/i7 c2 ) 1 / 2 , where H c2 is the upper critical field 
and N( 0) is the electronic DOS at the Fermi energy in 
the normal state. This pair breaking effect arises from 
Doppler shifts in the energies of ±k and — k states caused 
by the superfluid flow in the vortex state. For low H, 
pairs are broken near the nodes, where the s/c gap is 
small, but the region widens as H is increased. The effect 
is seen in heat capacity studies of YBa 2 Cu 307 crystals, 
for example Ref. H3J As shown in Fig. 0a), the low bias 
data show very little 17-dependence. This is presumably 
because of the strong 9 dependence for coherent c-axis 
tunnelling as explained in Refs. [T5J IS In contrast, the 
lower dip is strongly suppressed by fields of 13T applied 
perpendicular to the Cu0 2 planes and as shown later in 
Fig. 0 a), for OD80 it varies as A + BH 1 ? 2 , where A and 
B are fitting parameters. This must mean that the lower 
dip is strongly affected by the presence of pairs in the 
nodal regions. When they are broken by the applied field 
then the lower dip is attenuated. In contrast, as shown 
in Fig. 0b), dl/dV at low V varies approximately as 
A+BH . The 17-dependence of the upper dip is assessed 
by plotting dl/dV at 125 mV vs. H in Fig. 0c). Here the 
two OD mesas show remarkably similar behavior while 
the effect for OP86 is a factor of two larger. Differences 
between the two OD mesas at lower V in Figs. 0 a) and 
(b) are probably caused by the sensitivity of the nodes 
to inter- and intra-layer scattering processes. 

Fig. 0a) shows the overall T-dependence of our raw 
dl/dV data for OD80 in zero field at selected tempera- 
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FIG. 2: Color online: (a) dl/dV curves for OD80 at selected 
values of T, see the SM for the complete data set. (b) In¬ 
dependence of the d -wave gap 2Ao(T) up to T c from the main 
peaks in dl/dV for OD80 (black circles), OD78 (red squares) 
and OP86 (blue triangles). For OD80, green squares above 
T c =80 K give voltages of broad maxima in dl/dV. Green 
dashed lines show their increased breadth by marking regions 
where dl/dV > 0.95 (dl/dV) M Ax- 


tures. The data in Fig. [2^a) and for all three mesas in 
Fig-! show that the dip and the hump at higher V, are 
strongly T-dependent and have almost disappeared at 
50 K even though Aq has hardly changed from its low T 
value there. The attenuation of the hump is much smaller 
up to 40 K, but it shifts up with increasing T and also dis¬ 
appears rapidly between 50 and 60 K. It seems that this 
strong T-dependence, especially the shifts of the dips and 
humps with T, is unlikely to be caused by a conventional 
phonon pairing mechanism. As shown in Fig.[3j the lower 
dip is partially suppressed by a magnetic field, but it is 
not shifted, unlike the effect of temperature. For all three 
mesas a relatively sharp fall in 2Ao(T) also sets in just 
above 50 K as shown in Fig. [2jb). A striking feature of 
Fig. [2jb) is that just below T c , 2A 0 (T) ~ fi = 5.4fcsT c , 
possibly suggesting that the integrity of the magnetic 
mode is essential for superconductivity [T9] . 

The dl/dV curves in Fig. [2^a) above T c also have 
peaks whose breadth increases rapidly with T as indi¬ 
cated by the green dashed lines for OD80 in Fig. [2jb) . 
They appear to be states-conserving, for example at 84 K 
the polynomial normalization used earlier gives a dl/dV 
curve that conserves states to within 2.2% for the range 
of V shown in Fig. [2ja). The presence of these broad 
peaks for p = 0.19 agrees with the recent laser ARPES 
study of Bi-2212 showing a trisected s/c dome [8. and a 
pseudogap above T c extending up to p = 0.22. But at 
low T the pseudogap sets in abruptly below p = 0.19 [8] 
in agreement with earlier heat capacity [24. and penetra¬ 
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FIG. 3: Color online: Zoom of the structure above the gap 
edge, for (a) OD78, (b) OD80 and (c) OP86 at selected T. 
The data have been normalized for clarity, see Footnote 1T8I 
Full data sets including un-normalized data are given in the 
SM. Data at 1.4 K for a magnetic field of 13 T || c are also 
shown. 

tion depth [23 measurements. We note that the energies 
of the broad peaks above T c are comparable with 2Ao, 
so if they are caused by a competing mechanism [8] this 
would also tend suppress superconductivity once 2Ao(T) 
falls below a certain value. 

As recognized previously [9], the T-dependence of the 
data in Fig.[2](a) at low bias cannot be ascribed simply to 
thermal broadening. This is shown explicitly in Fig. |4jd) 
where the effect of thermal broadening has been calcu¬ 
lated by smoothing the dl/dV vs. V data at 1.4 K over 
a voltage window of 4 kgT/e. We therefore conclude that 
the larger values of dl/dV at low V, i.e. (dI/dV) res are 
caused by T-dependent pair-breaking processes in line 
with early break junction work [7]. The magnitude of the 
scattering rate (l/r p h) estimated from the energy broad¬ 
ening of the curves at low bias in Fig.|4](d) is surprisingly 
large, h/r p b ~ 14.5 meV or 170 K at 50 K. If real electron- 
boson scattering processes are responsible for 1 / r p b then 
its T-dependence will be related to the boson DOS. In 
Fig. [4](e) we plot (dI/dV) res vs. T. Normalized data for 
the three mesas are in excellent agreement, the data for 
OD78 and OD80 obey an A + BT 4 law, while for OP86, 
A + BT 3 gives a marginally better fit. Finally in Fig.^f) 
we show that the amplitude of the dip-hump feature for 
OD80 obeys an A — BT 2 law to high accuracy. The am¬ 
plitude for OD78 is smaller and is not such a good fit to 
A — BT 2 , probably because of non-ideal behavior near 
the nodes. 

In summary we have reported intrinsic SIS planar tun¬ 
nelling data for three crystals of the cuprate superconduc- 


























4 



FIG. 4: Color online: (a) normalized minima near 75 mV 
(lower dips) vs. H for the three mesas at 1.4 K, solid lines 
show fits to A + BH 1 ^ 2 , where A and B are fitting parame¬ 
ters. (b) dl/dV at low bias (10 mV) normalized to high bias 
(148 mV) vs. H. (c) dl/dV at a fixed voltage near upper 
dip (125 mV) vs. H , dashed lines are guides to the eye. (d) 
Zoom of some of the data for OD80 at low bias. Solid symbols 
show measured values, lines show the 1.4 K data smoothed by 
4.ksT to see the effect of thermal broadening. The horizon¬ 
tal dashed line shows that at 50 K the energy broadening is 
14.5 meV. (e) T dependence of the residual value of dl/dV at 
low V, which is caused by pair breaking, normalized to dl/dV 
at 148 mV for the three mesas. Fits to A + BT A (solid lines) 
and A + BT 3 (dashed line) are also shown, (f) T-dependence 
of the maximum in dl/dV at the “hump” minus the mini¬ 
mum at the “dip”, for OD80 and OD78, showing A — BT 2 
behavior, data for OP86 are less smooth and are not shown. 


tor Bi-2212, and discussed their field- and temperature- 
dependence. We argue that the dip structure must be 
connected with the nodes, presumably via the Eliashberg 
equations, because it is destroyed by H or T-induced 
pair-breaking. Accounting for this fact in terms of the 
magnetic mode that has Q ~ (7r/a,7r/a) [191 EH will be 
a strong test of this pairing model. T-dependent inelas¬ 
tic scattering is strong in the nodal regions, and could 
possibly be caused by the same excitations whose virtual 
exchange is providing the pairing “glue”. 
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FIG. 5: Supplemental Material. Typical I — V curves for the three mesas taken while sweeping the current up and down in a 
controlled manner, red points show data, blue lines are fits. Switching to another branch occurs when the critical (Josephson) 
current of a particular junction is exceeded and a voltage develops across it. The current is then swept down to a finite value 
before being increased again. Eventually, all Josephson currents are suppressed and the extreme right hand red curves, with the 
largest values of V are obtained. These are the same curves as those measured on downward sweeps, typically from V = 1.5 V 
to 0.05 V, and for which dl/dV is measured using a lock-in technique. The blue lines show fits of the form I = aV + bV 3 + cV s 
to the (N-l)th curve which are then scaled by n/(N-l) where n is an integer between 1 and N, for comparison with the 
measured values shown by red points. Differences between the red data points and the blue lines give an indication of possible 
non-uniformity in junction areas. 
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FIG. 6: Supplemental Material. This Figure provides confirmation of the statements at the bottom of the first column of page 
2 of our paper. The log-log plots show that the dl/dV data at 1.4 I< and 0 T for all three mesas obey the law A + BV ' 4 , 
below 20 mV or 0.74Ao. In the above figure the residual values A which have been subtracted are 0.309, 0.203 and 0.117 
x 10 ~ 3 Q~ 1 for OD78, OD80 and OP86 respectively. The calculated curve for completely incoherent tunnelling (no k 

conservation) between two weak-coupling d-wave superconductors is also shown. It has a completely different shape and a 
limiting V 2 dependence below V/N = 0.7 Aq. 
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FIG. 7: Supplemental Material. This Figure shows T-dependent tunnelling data for all three mesas at all values of T (the data 
for OD80 at 1.4 K are shown in Fig.2 of the paper). The areas of the mesas were measured by scanning electron microscopy. 
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FIG. 8: Supplemental Material. This Figure shows the tunnelling data at 1.4 K in fields from 0 to 13 T applied perpendicular 
to the Cu0 2 planes, for mesas OD78 and OP86 together with zooms of the lower and upper dips. The equivalent data for mesa 
OD80 are shown in Figure 1 of the paper. 
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FIG. 9: Supplemental Material. This figure shows the same data as in Fig. 3 of the paper but without normalization. 












